Astronomy & Astrophysics manuscript no. qwl 


© ESO 2008 


February 4, 2008 





Tracking quintessence by cosmic shear 

Constraints from VIRMOS-Descart and CFHTLS* and future prospects 

Carlo Schimd 1 ' 2 , Ismael Tereno 2 3 , Jean-Philippe Uzan 2 , Yannick Mellier 2 4 , 
Ludovic van Waerbeke 5 , Elisabetta Semboloni 2 , Henk Hoekstra 6 , Liping Fu 2 , and Alain Riazuelo 2 

VO "■ 

O ' DSM/DAPNIA, CEA/Saclay, 91 191 Gif sur Yvette cedex, France 

C> 2 Institut d'Astrophysique de Paris, UMR7095 CNRS, Universite Pierre & Marie Curie - Paris, 98 bis bd Arago, 75014 Paris, France 
3 Departamento de Fisica, Universidade de Lisboa, 1749-016 Lisboa, Portugal 
<-h ! 4 Observatoire de Paris - LERMA, 61 avenue de I'Observatoire, 75014 Paris, France 

Department of Physics and Astronomy, University of British Columbia, 6224 Agricultural Road, Vancouver V6T 1Z1, Canada 
Department of Physics and Astronomy, University of Victoria, Victoria V8P 52C, Canada 



6 



o 



I 

o 



Received February 4, 2008; accepted ... 

ABSTRACT 



> 

00 . 

I/-) . Context. Dark energy can be investigated in two complementary ways, by considering either general parameterizations or physically 
r- H ■ well-defined models. This article follows the second route and explores the observational constraints on quintessence models where the 
C^*) ' acceleration of our universe is driven by a slow-rolling scalar field. The analysis focuses on cosmic shear, combined with type la supernovae 
data and cosmic microwave background observations. 

Aims. This article examines how weak lensing surveys can constrain dark energy, how they complement supernovae data to lift some 
degeneracies and addresses some issues regarding the limitations due to the lack of knowledge concerning the non-linear regime. 
Methods. Using a Boltzmann code that includes quintessence models and the computation of weak lensing observables, we determine the 
O j 1 shear power spectrum and several two-point statistics. The non-linear regime is described by two different mappings. The likelihood analysis is 
based on a grid method. The data include the "gold set" of supernovae la, the WMAP-1 year data and the VIRMOS-Descart and CFHTLS-deep 
and -wide data for weak lensing. This is the first analysis of high-energy motivated dark energy models that uses weak lensing data. We explore 
larger angular scales, using a synthetic realization of the complete CFHTLS-wide survey as well as next space-based missions surveys. 
£-j Results. Two classes of cosmological parameters are discussed: () those accounting for quintessence affect mainly geometrical factors; ii) 
• • i cosmological parameters specifying the primordial universe strongly depend on the description of the non-linear regime. This dependence is 
, ^ addressed using wide surveys, by discarding the smaller angular scales to reduce the dependence on the non-linear regime. Special care is 
, payed to the comparison of these physical models with parameterizations of the equation of state. For a flat universe and a quintessence inverse 
% I ■ power law potential with slope a, we obtain a < 1 and Qqo = 0.75^^ at 95% confidence level, whereas a = 2+' 8 , £2qo = 0.74*qSj when 
including supergravity corrections. 
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1. Introduction 
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* Based on observations obtained with MegaPrime/MegaCam, a 
joint project of CFHT and CEA/DAPNIA, at the Canada-France- 
Hawaii Telescope (CFHT) which is operated by the National Research 
Council (NRC) of Canada, the Institut National des Sciences de 
l'Univers of the Centre National de la Recherche Scientifique (CNRS) 
of France, and the University of Hawaii. This work is based in part on 
data products produced at TERAPIX and the Canadian Astronomy 
Data Centre as part of the Canada-France-Hawaii Telescope Legacy 
Survey, a collaborative project of NRC and CNRS. 



Cosmological observations provide increasing compelling ev- 
idences that the expansion of the universe is accelerating 
and that the cosmic history of the universe seems today 
dominated by another component than its matter and radi- 
ation content (see e.g. IPeebles & Ratra 20031 ICarroll 200T1 
Padmanabhan 20031 IPeter & Uzan 20031 chap. 12 for reviews 
and references therein). If so, one of the most challenging issue 
of fundamental physics is to understand the cause of this accel- 
eration, a question often referred to as the nature of the dark 
energy. Various solutions, from the introduction of a new type 
of matter to a modification of general relativity to describe the 
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gravitation interaction, have been considered. A classification 
of these models with some relevant observational tests that can 
help to distinguish between each class from their underlying 
new physics is discussed in Uzan (2004) and Uzan, Aghanim 
& Mellier (2004). 

Dark energy appears in the Friedmann equations through 
its effective density and pressure. Data are usually interpreted 
assuming the validity of the Copernician principle (so that the 
dynamics of spacetime is completely described by a single 
function, the scale factor a) and the validity of Einstein equa- 
tions (and thus the standard Friedmann equations), so that the 
density and pressure of the dark energy component are de- 
fined by p de = (3/SkG)(H 2 + K/a 2 ) - p m - p r and P de = 
(-l/8nG)(a/a + H 2 + K/a 2 ), where H = a/a is the Hubble 
parameter and a dot refers to a derivative with respect to the 
cosmic time, while p m and p r are the density of pressureless 
matter and radiation, respectively. It follows that the equation 
of state of the dark energy corresponds to 

3D. de w = -1 +fl f + 2q, (1) 

where q = -aa/a 2 is the deceleration parameter and = 
-K/a 2 H 2 . From this point of view, w characterizes the dynam- 
ics of the cosmic expansion. More precisely, it parameterizes 
the deviation, H(z) - H(z), between the Hubble function of 
the observed Universe, H(z), and that predicted for a universe 
filled only with pressureless matter and radiation, H(z). It is 
therefore equivalent to specify w(z) or H(z) - H(z). However, 
when general relativity is assumed to describe gravity, w re- 
duces to Pde/Pde so that, in addition to the deviation from H(z), 
it also gives some insight on the properties of dark energy (see 
e.g. Martin, Schimd & Uzan (2005) for a case in which w does 
not reduce t the equation of state of a matter component). 

Although effective equation of state derived from obser- 
vations is a key empirical information on the rough nature of 
dark energy, a detailed description of its properties demands 
more thoughtful data interpretation. For example, all geometri- 
cal observables rely on the integration of the Hubble parameter, 
hence on a double integration of the equation of state w, that 
eventually dilutes or totally washes out its possible redshift de- 
pendence. If w is close to -1, as observations tend to indicate, 
then it is in general difficult to demonstrate by geometrical tests 
that w + -1 or that dw/dz + 0; both would exclude a pure cos- 
mological constant. Exploring early properties of dark energy 
models would be even more challenging since for w - — 1 the 
ratio between the matter and dark energy densities scales ap- 
proximatively as (1 + z) 3 so that dark energy is dynamically 
negligible at redshift z > 2. It leaves little freedom to deter- 
mine the scaling of the dark energy density and to demonstrate 
that it is not properly described by a power law [(1 + z)"], as 
would be the case for a constant w (see e.g. |Kujat et al. 2002) . 

From the theoretical point of view two routes can be fol- 
lowed. One can either exhibit a general "model-independent" 
parameterization of the equation of state of the dark energy, as 
discussed in the previous paragraph, or rely on a completely 
specified theoretical models. A useful parameterization has to 
be realistic, in the sense that it should reproduce predictions 
of a large class of models, it has to minimize the number of 
free parameters and to be simply related to the underlying 



physics (see e.g. Linder & Huterer 2005). Because the result of 
the analysis will necessarily have some amount of parameteri- 
zation dependence (Basset et al. 2002 1, choosing the specified 
physical model strategy seems preferable to break degenera- 
cies. In particular, it enables to compute without any ambiguity 
their signature both in low and high redshift surveys, such as 
the cosmic microwave background (CMB). The increasingly 
flourishing number of models hampers to provide a compre- 
hensive set of unambiguous predictions to constrain physical 
models one by one with present-day observations, but there are 
still several benefits in exploring dark energy this way, in par- 
ticular when weak lensing surveys are used together with CMB 
observations. This is deeply related to the evolution of dark en- 
ergy properties and the growth rate of structure with look-back 
time, as discussed below. 

At low redshift, w suffices to get observable that are 
all functions of H(z) (see for example Peebles 1993, 
IPeter & Uzan 20051 . This is the case of all background 
quantities (e.g. luminosity distance, angular distance, look- 
back time, etc.) as well as of the linear growth factor of 
density perturbations. It follows that the equation of state 
encodes all relevant information, provided the amplitude 
of the power spectrum is calibrated by adding a new pa- 
rameter, cr 8 , the variance of the density perturbation on a 
scale of 8 ft -1 Mpc (see e.g. IBenabed & Bernard eau 20011 
IBenabed & van Waerbeke 20031 IDoran et al. 20001 

However, as far as weak lensing is concerned, 
it was shown JBenabed & van Waerbeke 2003 ; 
IBenabed & Bernardeau 200lT > that, for a fixed redshift of 
the sources, the modification of the growth factor in the linear 
regime was degenerate with the normalization factor. Hence 
normalizing on the CMB avoids this problem, and at the same 
time it is important to describe the non-linear regime. The 
use of CMB together with weak lensing data is therefore a 
logical way to constrain specified theoretical models, beyond 
the description by an empirical equation of state. 

At higher redshift, and in particular to relate the amplitude 
of the matter density power spectrum to the one of the primor- 
dial power spectrum, one would need to include a description 
of the evolution of the perturbations of dark energy. In particu- 
lar, this effect becomes increasingly important as w approaches 
zero ( Benabed & Bernardeau 2001 ). Note that dark energy per- 
turbations have a non-adiabatic component that also requires a 
detailed model to be described. This depends on the physical 
model of dark energy and cannot be incorporated in a simple 
model-independent way. 

While the ability of lensing data to constrain the 
equation of state of dark energy has been widely stud- 
ied llBenabed & van Waerbeke 20031 IHu & Jain 20041 
|Jain & Taylor 2003| >, there have been very few analysis 
with real data. Hoekstra et al. (2005) and Semboloni et al. 
(2005) used the CFHTLS wide and deep data to constrain 
a constant equation of state. Jarvis et al. (2005) analyzed 
the 75 square degrees CTIO lensing survey, combined with 
type la supernovae (Sn la) data and CMB, assuming a constant 
equation of state and a parameterization of the form proposed 
by |Chevallier & Polarski (200T)| and then by |Linder (20031 
In this article, we consider a class of completely defined 
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quintessence models, realized by a self-interacting scalar field. 
Hence all observational signatures (Sn la, lensing, CMB) 
can be explicitly computed, with no ambiguity in the way to 
deal with the perturbations of dark energy. This theoretical 
extension of the standard ACDM model involves only one ad- 
ditional parameter, needed to characterize the self-interacting 
potential of the quintessence field. It follows that, as discussed 
above, we will be able to normalize our initial power spectrum 
on the CMB angular power spectrum and, as a consequence, 
erg and any possible dependence on the shape of the analytical 
fit of the transfer function will disappear from our discussion; 
the value of <x g will be an output of each models. The problem 
of the pivot redshift (see § 12. 2> that appears when combining 
different datasets also disappear in that approach. As a con- 
clusion, this approach is very efficient in terms of the number 
of extra-parameters and of the interpretation of the data. Let 
us emphasize that, even though we also consider Sn la and 
CMB data, we will focus on weak lensing - cosmic shear data. 
This article, being the first analysis of the CFHTLS data for 
dark energy studies, illustrates the power and the problems of 
lensing survey in studying dark energy. 

The article is organized as follows. In Section[2] we define 
the quintessence models we are considering and we recall their 
main properties. We also compare them to various parameter- 
izations proposed in the literature. Section [3] focuses on cos- 
mic shear. After a reminder on theoretical issues, we describe 
the weak lensing data used for our analysis and the computa- 
tional pipeline, finally we outline the likelihood analysis on real 
(§ 13. 5> and synthetic (§ 13. 6i datasets. In Section|4]we combine 
weak lensing data with type la supernovae and CMB temper- 
ature anisotropies; Figure^] and Table ^ summarize the con- 
straints on quintessence parameters. We finish, in Section|5] by 
an estimation of the proficiency of two possible space-based 
wide field imagers to unveil the nature of dark energy. 



2. Modelling dark energy 

2.1. Quintessence models 

In this work, we consider the simplest class of quintessence 
models (Ratr a & Peebles 19881 IWetterich 1988i in which a 
scalar field Q is slow-rolling in a runaway potential. Numerous 
forms of potentials have been proposed but we restrict to two 
classes of potentials. Let us briefly summarize the properties of 
these models. 

The first class of potentials is an inverse power law 
JRatra & Peebles 19881 IWetterich 19 88 1 



V(Q) = M 4 (Q/M p y 



(2) 



often quoted as Ratra-Peebles (RP) potential. M p = (8nGy 1/2 
is the reduced Planck mass. The potential depends on two free 
parameters: a is a positive index while M is a mass scale that 
has to be adjusted to fit Oq today, once a is given. In particular, 
if Qq dominates today then M and a are related by 



log 



(i 



M \ 
GeW 



19a - 47 
4(a + 1) ' 



(3) 



The second class of potential is an extension of the previ- 
ous that takes supergravity (SUGRA) corrections into account 
when Q ~ M p dBrax & Martin 1999i . 



V(Q) = M 4 (£/M p )- a exp(£ 2 /2M p 2 ). 



(4) 



Both potentials have a similar dynamics as long as Q <K M p but 
differs at low redshift, in particular concerning their equation 
of state. In the SUGRA case, it is pushed toward -1 and one 
expects w ~ -0.82 (jBrax & Marti n 1999i . 

With such well-defined models, the dynamics of the back- 
ground is completely characterized by the Klein-Gordon equa- 
tion for the scalar field, 



dV 

Q + 3HQ + — = 0, 
dQ 

in addition to the Friedmann equation 



(* 2 +S) 



8kG 



Q 2 



(5) 



(6) 



allowing for this new matter contribution. These equations 
characterize background and low redshift observations and in 
particular the linear growth factor. Let us stress that, at this level 
of description, one can describe the quintessence component as 
a fluid with a time-dependent equation of state. This is due to 
the fact that the speed of sound, c 2 , is given by 



2-1 i-L^L 

° s 3 HQ dQ 

and that the equation of state evolves as 
w = -3H(1 + w)(c] - w). 



(7) 



(8) 



These models share the interesting property to possess scal- 
ing solutions which are attractor of the dynamical evolution. In 
general, but depending on the initial conditions, the dynamics 
starts with an early kinetic phase (Q 2 » V) in which w ~ 1 so 
thatpQ oc (1 +Z) 6 . The field behaves as Q = Q, < -A{\ +z) and it 
freezes to a constant value. Since the kinetic energy decreases 
while the potential remains constant, this regime cannot last 
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Fig. 1. Evolution of the equation of state w and of the sound 
speed c s with the redshift for an inverse power law quintessence 
model with a — 6, including or not the supergravity correction. 
We recover, from high to low redshift, the kinetic, slow-rolling 
and tracking phases described in the text. 
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forever. When the potential starts dominating, the equation of 

state shifts suddenly to w 1 . During this transition regime, 

pq oc (1 + z) . Then, there is a potential regime that lasts until 
the tracking regime during which 



a + 2 



(9) 



where wb stands for the equation of state of the fluid domi- 
nating the background. At that stage, the scalar field is slow- 
rolling so that w < 0. Fig.^depicts the evolution of w and c 2 
during these various regimes. 

When cosmic microwave background anisotropics and 
large scale structures are considered, one needs to include the 
description of the evolution of the perturbations, and in partic- 
ular include those of the scalar field described by 



. (k 2 d 2 V\ 
6Q + 3H6Q + l [ - 2+ — 2 )60-S = i) 



(10) 



where S encodes the perturbations of the metric of the space- 
time and k is the comoving wavenumber of the perturbation. 
It was shown (IBrax et al. 20001 IRiazuelo & Uzan 20021 that 
there exists an attraction mechanism for super-Hubble wave- 
length so that the spectrum is insensitive to the initial condi- 
tions for the scalar field. 

2.2. Models and parameterizations 

Most data, and in particular supernovae data, are being an- 
alyzed using a general parameterization of the equation of 
state. These parameterizations are useful to extract model- 
independent information from the observations but the inter- 
pretation of these parameters is not always straightforward. In 
this paragraph, we remind the properties of some interesting 
parameterizations and compare them to the quintessence mod- 
els we are considering. 

Let us first recall that general parameterizations of the equa- 
tion of state as 



w(a) = w(ao) + [w(a m ) - w(«o)] T(a, a,, A) 



(11) 



were shown to allow an adequate treatment of a large 
class of quintessence models ( |Corasaniti & Copeland 2 003 
Basset et al. 2002 1. Such a parameterization involves four pa- 
rameters {w(flo), w(a m ), a t , A} and a free function F varying 
smoothly between one at high redshift to zero today. Even 
though it reproduces the equation of state of most quintessence 
models, it is not economical in terms of number of parame- 
ters since most quintessence potentials involve one or two free 
parameters. If one assumes that the parameterization is sup- 
posed to describe the dynamics of a minimally coupled scalar 
field, the knowledge of w is sufficient but in a more general case 
one would need more information: The background dynamics 
depends on the potential and its first derivative, which can be 
related to w and w. Accounting for perturbations, one needs to 
know the second derivative of the potential [see Eq. (II 0H which 
can be inferred from w (Dav e et al. 2002> . 

Since we expect dark energy to have observable conse- 
quences on the dynamics only at late time, one can consider 
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Fig. 2. Deviation of quintessence equation of state for Ratra- 
Peebles (solid) and SUGRA (dotted) models with a — 6 from 
the generalized parameterization, Eq. (II 1> . setting z p i vot = 
(thick) or z p i vot = 0.5 (thin). Fitting the previous one up to 
z < 0.3, a deviation larger than 2% occurs at z - 1 for Ratra- 
Peebles models while at z — 0.5 for SUGRA models. 



an equation of state obtained as a Taylor expansion around a 
pivot point, 



w(a) = Wp ivot + w a (a pivot - a). 



(12) 



This form depends on only two parameters and is 
a generalization of the parameterization proposed 
by |Chevallier & Polarski (200T)| and then |Linder (20031 
where a p i vo t = 1 • Two considerations are in order when using 
such a parameterization. First, the redshift band on which this 
is a good approximation of the equation of state is unknown. 
Clearly, compared with the form (II Q , it is unlikely to describe 
dark energy up to recombination time; see Fig. |2] Secondly, 
when combining observables at different redshift such as weak 
lensing, Sn la and CMB, one should choose the value of a p j vo t 
in such a way that the errors in w p i vo t and w a are uncorrected 
(IHu & Jain 2004> . It follows that the pivot redshift is the red- 
shift at which w is best constrained. In particular it was argued 
that it is important to choose a p i vo t + 1 for distance-based 
measurements. The problem lies in the fact that the pivot 
redshift is specific to the observable. In this respect, dark 
energy models defined by a Lagrangian are more suitable, 
yielding to a definite equation of state as a function of redshift, 
hence more general than a Taylor expansion around a pivot 
point. Eventually, one can read out the values of w p j vo t and w a 
at whatever redshift. Fig. [3] depicts the value of w p i vo t and w a 
for the quintessence models we use, Eqs. (|2j and (|4j, for the 
pivot redshifts z p i vot = and z p i vot = 0.5. This complication, 
arising when one wants to combine datasets with different 



Zpivot- 



will also make it more difficult to infer constraints on the 



physical models from the constraints on the parameterization. 

There is an alternative way to get a first hint on the nature 
of dark energy. It may be useful to consider the plane (w, w') 
where w' = dw/dlna is the derivative of w with respect to the 
number of e-folds. It was recently shown by Caldwell & Linder 
(2005) and Scherrer (2005) that quintessence models occupy a 
narrow part of this plane. This can be understood from Eq. (|S} 
which implies that w' + 3(1 - w 2 ) = 3(1 + w)(l - c 2 ). For 
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Fig. 3. Contour plots of the quintessence equation of state. We compare the equation of state of two quintessence models with 
the parameterization d!2> for two values of the pivot redshift: z p i vot = (left) and z p i vot = 0.5 (right). Solid lines correspond to 
level contours for w p j V ot while dotted lines correspond to level contours of w a . We have chosen the spacing of all the contour lines 
to be Aw = 0.1, except for the plots in the upper line, where Aw a = 0.02. The upper line corresponds to Ratra-Peebles models, 
Eq. 0, while the lower line corresponds to SUGRA models, Eq. 0). Due to the exponential correction, wq is always smaller for 
SUGRA models because the potential is flatter and the field is rolling slower. Also, the value of w p i vo t and w a are more sensitive 
to the choice of z p j vot for SUGRA models than for Ratra-Peebles models. 



quintessence models, 1 + w > 0, and Eq. @ implies that c] < 1 
(because Q > and V < 0) so that 



w > -3(1 



■ w 2 ), 



(13) 



without any assumptions on the dynamics of the scalar field. In 
Caldwell & Linder (2005), two classes of quintessence models 

where exhibited, namely "thawing" models, in which w 1 

initially and increases as Q rolls down the potential , and 
"freezing" models, in which w > — 1 initially and tends to- 
ward - 1 as Q rolls down the potential. "Freezing" models con- 
tain tracking models and in particular the Ratra-Peebles mod- 
els, Eq. 0, and SUGRA models, Eq. ®, considered in this 
work. Using a combination of numerical simulations and phys- 
ical arguments, they concluded that 



3w(l + w) < w' < 0.2w(l + w) 



(14) 



for "freezing" models. From an observational point of view, 
the analysis of supernovae data (Riess et al. 2004) showed that 
if the universe is flat then w = -1.31+° ?J and w' Q = -1.48+J}* 
after marginalizing on Q. m Q. If one further imposes that wq > - 1 
then wo < -0.76 and w' = -0.6 + 0.5. In the case where w is 
assumed constant then w = -1.02*SHq and w < -0.72 at 68% 
and 99% confidence level respectively. 

This phase-space analysis creates a link between the var- 
ious parameterizations and the physical models. It can be 
shown JScherrer 20 05 1 that different classes of models (e.g. 
k-essence, Chaplygin gaz, quintessence, etc.) lies in different 
parts, hence offering a way to distinguish between these models 




Fig. 4. Dynamics of the two quintessence models in the plane 
(w, w'). The shaded regions correspond to the constraints Jl 31 
in light gray and J 141 in dark grey. We have considered a 
Ratra-Peebles (solid) and SUGRA (dash) models with a = 6 
(thick/red) and a — 1 1 (thin/blue). Only in the tracking regime 
the models are compatible with Eq. dl4t . 



without measuring w{z). In particular, Fig.@]depicts the dynam- 
ics of some Ratra-Peebles and SUGRA models in the (w, w') 
plane, superposed to the regions where the inequalities Jl 3i 
and (TBI hold. Notice that the trajectories for the SUGRA po- 
tential are essentially the same of the Ratra-Peebles one, de- 
viating just at low redshift towards the cosmological constant 
solution. 
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2.3. Summary 

Quintessence models require only two parameters to describe 
the whole dynamics (with no redshift limitation). Compared to 
a pure cosmological constant, described by only one number, 
this gives us one extra parameter. In terms of extra-parameter 
with respect to a standard ACDM, this is equivalent as consid- 
ering a constant equation of state for dark energy. 

On the other hand, a parameterization of the equation of 
state of the dark energy is sufficient to describe the low redshift 
universe. But, the parameterizations which describe accurately 
quintessence models involve at least four extra-parameters. It is 
thus more economical to work directly with the physical model. 
In that case, the evolution of perturbations can be inferred from 
w but this in not the case in more general situations. 

Parameterizations with fewer parameters have a smaller dy- 
namical range and are likely to be bad descriptions of dark en- 
ergy at high redshift, notably for CMB computation. Besides, 
when combining dataset, the choice of the pivot redshift may 
induce some difficulties in interpreting the constraints on the 
parameters in terms of constraints on the physical models. 

So considering directly a well-defined physical model in- 
stead of a parameterization is economical concerning the num- 
ber of extra-parameter and avoids the problem of the pivot red- 
shift. It allows us to compute the prediction of the models at all 
redshift (Sn la, weak lensing, and CMB). On the other hand, it 
concerns only a small class of models. 

This discussion shows that both routes are complementary. 
In particular, it would be worth to evaluate to which accuracy 
constraints of order 1 % on a given parameterization constrain 
physical models. 



3. Cosmic shear 

Gravitational lensing by large scale structures of 
the universe produce weak distortion fields and col- 
lectively modify the shape of background galaxies 
(see e.g. IBartelmann & Schneider 20011 IMellier 19991 
jRefregier 2003) . Though this is a very weak signal and 
a challenging task, it has been detected almost simulta- 
neously by van Waerbeke et al. (2001), Wittman et al. 
(2000), Kaiser et al. (2000), Bacon et al. (2000) and it now 
routely observed by many groups around the world (see 
lyanW aerbeke & Mellier 20031 IHoek stra 2003 for recent 
reviews on observations). Over the past years the huge efforts 
carried out by these groups to deal with critical systematics 
considerably improved the reliablity of the lensing signal 
and stenghtened the ability of cosmic shear to constrain 
cosmological models from the statistical analyses of galaxy 
ellipticies. 



3.1. Generalities 

The gravitational lensing effect depends on the second order 
derivatives of the gravitational potential. The convergence, k, 
and the shear, y = (71,72), describe the distortion of back- 



ground images by the matter along the line of sight. These com- 
ponents are related by 



Ak 



(d? - 5^)71 + 23^272. 



(15) 



The evolution of the convergence is dictated by the Sachs 
equation ( Sachs 1962 see lUzan & Bernardeau 200"0l for a more 
modern description). The convergence in the direction 6 can be 
related to the matter distribution integrated along the line of 
sight 



3 £2m0 

2 



Sk(x -x')Sk(x') 
Sk(x) 

j m [s K (x')e, x r 

a(x') 



dx' 



for sources located at a radial distance x defined by 
dz' 



X(z) 



-f 

Jo 



H(z') 

and S k is the angular diameter distance given by 

sin( \fK X )/ <K K>0 

S K (X) = \x K = . 

sinh( ^PKx)I^-K K < 



(16) 



(17) 



(18) 



If the sources have a distribution given by n x (x)dx = «(z)dz 
then the effective convergence takes the form 



<€) = 



Xll 



nJx)K{0,x)dx 



(19) 



where xh is the comoving radial distance of the hori- 
zon. Decomposing the convergence in 2-dimensional Fourier 
modes, 



r d 2 e 
m - J -* 



(20) 



the shear power spectrum, defined by {k(t)k(t')) - P K ({)6 (2) ({+ 
('), can be related to the 3-dimensional power spectrum of mat- 
ter density perturbations P m by 



9 H*Q? n 



f 

Jo 



\g(xY 


2 


\S K (X) 1 


a(x)_ 




[ I * 



dx 



(21) 



in the small angle approximation (see e.g. 
IBartelmann & Schneider 20011 IPeter & Uzan 20051 chap. 7). 
The function g is given by 



, v r m , „ sk(x'-x) a , 



(22) 



Note that the window function W(z) = [g(x)/ a (x)] 2 is peaked 
around z — z* 12 for a distribution of sources n(z) approxi- 
mately peaked at redshift z* ■ This will be useful in the choice of 
the pivot redshift. Let us stress that expressions d!6i and \2\\ 
assume the validity of the Poisson equation and thus of gen- 
eral relativity. These expressions may be slightly different in 
more general contexts and even be used to test general relativ- 
ity ( see lUzan & Bernardeau 2001l lSchimd et al. 20051 1. 

Neither k nor P K are directly observable, but only filtered 
quantities can be obtained. Cosmic shear can be measured by 
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various types of 2-point statistics which differ only by the cho- 
sen filtering scheme. This implies that their sensitivity to the 
power spectrum, and also to systematics, are different. In this 
work we consider two of these statistics, namely the aperture 
map variance, defined by 



{MM 



288 r 
U = — J «V 



(0 



J 4 (€6 C ) 



d{, 



(23) 



which is a bandpass estimate of the convergence power spec- 
trum, and the top-hat shear variance 



(r)(0c) 



d(, 



(24) 



which is a lowpass estimate of P K . Here J„ are the Bessel func- 
tions of the first kind. Both statistics can be deduced from 
two linear combinations of the radial and tangential com- 
ponents of the shear variance, £ ± = (y t 2 ) ± (y 2 ) (see e.g. 
Bartelma nn & Sch neider 2001 1, which are directly estimated 
from the shapes of background galaxies. 

3.2. Matter power spectrum and non-linear regime 

In the previous expressions, and in particular in Eqs. d!6i 
and (I21> . P m refers to the full 3-dimensional power spectrum 
of pressureless matter, including cold dark matter and baryons. 

In the linear regime, the growth factor D+ is the growing 
solution of 



D + + 2HD + - -H 2 Q. m D + = 0. 



(25) 



The second term that describes the damping due to the cos- 
mological expansion contains all the effect of dark energy 
on D + . In this regime, the total effect on lensing, Eq. \2\\ . 
is similar to the one obtained from a single redshift plane 
(see Benabe d & van Waerbeke 20031 ISchimd et al. 2005> so 
that the integrated growth effect is degenerated with the nor- 
malization of the spectrum. 

Dealing with low-redshift sources, the lensing pre- 
dictions and analysis involve the non-linear power 
spectrum. This regime cannot be described analytically 
from a perturbation approach (see however the proposal 
bv ICrocce & Scoccimarro 20 05 1 and one would need to rely 
on iV-body simulations. 

N-body simulations including quintessence were recently 
performed using the GADGET code (Dolag et al. 2003 1 or an 
adaptive refinement tree code ( |Klypin e7aT^003| . It was ar- 
gued by Dolag et al. (2003) that the halo concentration dis- 
tribution around the mean value does not depend on the cos- 
mology, while the concentration parameter depends on the 
dark energy equation of state at the cluster formation red- 
shift through the linear growth factor. However a systematic 
study confirming this claim is lacking. Klipin et al. (2003) and 
Mai nini et al. ( 2003) show that dark energy changes the virial 
density contrast, A c , which induces a change in the power spec- 
trum at small scales (k ^ 1/iMpc -1 ) but, for constant w, the 
error is smaller than the error in the expected non-linear model 
dJarvis et al. 2005> . 



Instead of specific A^-body simulations accomodating 
quintessence, in order to compute the non-linear matter 
power spectrum one can deal with linear-to-non-linear map- 
pings, for instance based on the stable clustering ansatz 
( Hamilton et al. 19911 IPeacock & Dodds 19961 or to a halo 
model (e.g. |Ma & Fry 20"00| |Seljak 2000| ISmith et al. 20031 . 
These mappings have been tested for several cosmologies in- 
cluding ACDM but not for dynamical dark energy models. 
Given their robustness and the precision level we can actu- 
ally reach, we can hope that they remain valid for the class 
of models we are considering. Indeed this is a very strong as- 
sumption that can be justified by the fact that we do not ex- 
pect the scalar field to cluster on small scales so that it is un- 
likely to affect the small scale behaviour of matter but by its 
influence on the expansion rate. Arguing that the clustering 
scale of the quintessence field is given by its Compton wave- 
length, Ma et al. (1999) propose an analytic approximation for 
the Peacock & Dodds (1996) formula to include quintessence 
with constant equation of state, claiming a 10% level accuracy. 
McDonald et al. (2005) propose a recipe to extend the afore- 
mentioned mappings to w + -1 cases, attaining a better ac- 
curacy for k % 10 Mpc~' by systematically exploring a wider 
parameter space. Nevertheless, upcoming weak lensing mea- 
surements require an improved description of the smaller scales 
physics (Huterer & Takada 2005 1, eventually including hydro- 
dynamics dZhan & Knox 20041 |jm"g et al. 2005) . 

In conclusion, it seems early to decide how dark energy, 
and quintessence in particular, modify the mapping calibrated 
on ACDM. For that reason, in this work we will consider two 
linear-to-non-linear mappings, by Peacock & Dodds (1996) 
and Smith et al. (2003), and try to identify the parameters that 
are not sensitive to this choice. We will also try to quantify how 
the other parameters are affected so that we can estimate how 
our ignorance of the non-linear regime limits the use of weak 
lensing. 

Hopefully, as we will also show, weak lensing data in the 
linear regime shall be able to be used. In that case, we can get 
an interesting constraint on the growth factor without messing 
with non-linear physics. 



3.3. Lensing data 

We use three sets of data for weak lensing, the VIRMOS- 
Descart ( Ivan W aerbeke et al. 2005t . the deep field of the 
CFHTLS survey ( Sem boloni et al. 20051 1 and the wide field of 
CFHTLS (IHoekstra et al. 20051 

The details regarding the survey properties, image and 
catalogue processing of VIRMOS-Descart data are described 
in Ivan Waerbeke et al. 200T1 and IMcCracken et al. 20031 The 
shear measurement and error analysis are described in 
Ivan Waerbeke et al. 200 5 and summarized in Table 2 of that pa- 
per. This survey covers an effective (the unmasked area) area of 
8.5 deg 2 spread in four fields. It probes lensed galaxies down to 
the limiting magnitude Iab = 24.5 has an effective galaxy num- 
ber density (after all selection processes) of 15 gal/arcmin 2 , and 
explores angular scales up to 50 arc-minutes. 
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CFHTLS 
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Fig. 5. Pipeline implemented for this work. Presently, we restrict to three free cosmological parameters, {Qqo, a, n s ), keeping 
fixed the others. The lensing code manipulates both background's and perturbations quantities computed by the CMB code, 
using a sources distribution n(z) depending on the used dataset. In particular, the source redshift parameter z s is left to vary and 
marginalized over afterwards. Finally, the likelihood is computed using (either real or synthetic) cosmic shear and Sn la data, 
both separately and jointly. The temperature CMB data are used to fix the amplitude A of the power spectrum at decoupling, and 
to put (conservative) constraints on the (£2qo, a) parameters sub-space using the location of the first peak. In the CMB section 
of the pipeline, we indicate by 5 m and <Sq energy density fluctuations in matter and quintessence components, respectively, while 
<$>(k, z) and *?(£, z) are the scalar perturbations of the metric (Bardeen potentials) in Fourier space. In the lensing section, we 
denote by me scale at which the power spectrum becomes non-linear, Pik^, z) ~ 1 . See §|5]for details. 



For the CFHTLS deep and wide, all relevant details are 
given in ISemboloni et al. 20051 and IHoekstra et al. 20051 re- 
spectively. The deep covers an effective area of 2.2 deg 2 in 
three fields, down to Iab = 26, but only samples angular scales 
up to 30 arc-minutes. The effective galaxy number density is 
22 gal/arcmin 2 . In contrast, the wide has a much larger effec- 
tive area than the deep (22 deg 2 ) but only spread in two fields. 
It explores angular scales up to one degree at about the same 
depth as VIRMOS-Descart. The effective galaxy number den- 
sity is 13 gal/arcmin 2 (the final selection produced a catalogue 
somewhat less deep than VIRMOS-Descart). Note that once 
completed the wide survey will be composed of three compact 
areas that will sample angular scales up to 5 degrees in three 
independent fields so that linear scales will be explored with 
much more accuracy than present-day CFHTLS wide data. 



Each survey used complementary photometric or spectro- 
scopic galaxy samples to derive the redshift distribution of the 
lensed galaxy samples. As discussed in the VIRMOS-Descart 



and CFHTLS cosmic shear papers, it is convenient to describe 
the redshift distribution by a three-parameter function 



n(z) = 



P 



-, exp 



(26) 



a and /3 are obtained from a fit of photometric redshift distri- 
butions derived from external redshift calibration surveys. The 
CFHTLS and VIRMOS are different lensing data sets, but they 
were obtained with the same telescope (although with differ- 
ent instruments) and with the same exposure time. Therefore 
the two lensing surveys have to be calibrated using the same 
external redshift data set in order to preserve the homogene- 
ity of the analysis. In van Waerbeke et al. (2005), VIRMOS 
was calibrated using the Hubble Deep Fields and MS 1008 (see 
van Waerbeke et al. 2001 for the details), while the CFHTLS 
lensing data (Hoekstra et al. 2006; Semboloni et al. 2006) were 
calibrated using the HDF only. The HDF and MS 1008 data in 
visible and near infrared bands produced accurate photomet- 
ric redshifts (Yahata et al. 2000; Athreya et al. 2002). We de- 
cided to use the HDF only for both VIRMOS and CFHTLS 
and abandon the MS 1008 field. We have checked that the re- 
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suits discussed in this paper do not depend whether or not we 
include the MS 1008 calibration field, since the sample variance 
due to the use of combined redshift calibration sets is absorbed 
in the redshift error (van Waerbeke et al. 2006). 

In the following, we will use for the CFHTLS-deep survey 
(Sem boionlet al. 2005> 



a = 1.9833, /? = 0.6651, z s 



n HQ8 1 +0.0129+0.0209 n7 x 
U,IWSI -0.01 14-0.0161 



giving a mean redshift (z) = 1.01, while for the CFHTLS-wide 
and the VIRMOS-Descart surveys dHoekstra et al. 2005> 



1.35, £= 1.654, z s 



n A/;s+0.035+0.053 
U - OUO -0.036-0.055 



(28) 



giving a mean redshift (z) = 0.76. We quote the lcr and 2<x er- 
rors. We use the same parameters for the VIRMOS-Descart and 
CFHTLS-wide redshift distributions since both surveys have a 
similar depth and the same effective galaxy number density. 

3.4. Description of the pipeline and choice of the free 
parameters 

The pipeline we developed is summarized in Fig. EI W e com- 
pute the evolution of background and perturbations power 
spectra in linear regime by means of a Boltzmann code al- 
lowing for photons, neutrinos, baryons, cold dark matter and 
quintessence scalar field. Notice that this code deals with 
several gauge choices and can also account for scalar-tensor 
theories of gravity (Riazuelo & Uzan 2002 1, hence allowing 
to study deviations from general relativity in this frame- 
work as well as extended quintessence scenarios dUzan 1 999 
lAmendola 2000llChiba 19 99 1. 

Using this code we compute the CMB temperature (TT) 
angular power spectrum (Q) in order to fix the amplitude of 
the initial matter power spectrum at the redshift of the last scat- 
tering. We do it by matching the computed Q with WMAP- 
lyr data (Hinshaw et al. 2003) at a high multipole, to be pre- 
ferred when studying quintessence since at low multipoles the 
integrated Sachs-Wolfe effect is dominant so that temperature 
anisotropies are not directly related to the density perturba- 
tions. Definitely, we use the multipole i ^ 110 of WMAP-lyr 
data, where the total relative error of the TT spectrum is smaller 
than 3%. This normalization procedure holds until the correla- 
tion between multipoles is weak and does not take into account 
the measurement errors on the amplitude of the TT spectrum. 
However, we expect that the final results would not be strongly 
affected by a more accurate normalization. Notice that, as men- 
tioned earlier on, <x 8 data is not used to normalize the spectrum. 
Hence, its value may be evaluatd from the matter power spec- 
trum and compared with the observed values. 

Once the linear matter power spectrum is known at ev- 
ery redshift, the weak lensing add-on code ( Schi md et al. 2 005 ) 
computes the non-linear power spectrum using two mappings, 
by Peacock & Dodds (1996) and Smith et al. (2003). We con- 
sider wavevectors ranging up to 10/; Mpc~' (corresponding to 
1 arcmin at z ~ 1 for cosmologies close to ACDM). By Limber 
projection, we deduce the shear power spectrum allowing for a 
source redshift distribution of the form i26\ . To finish, several 
two-point statistics in real space are computed, namely the top- 
hat shear variance (y 2 ), the aperture mass variance {M 2 ), as 



well as the two-point correlations £ ± . Let us stress that, like the 
Boltzmann code, the lensing add-on code works with scalar- 
tensor theories of gravity as well. 

Cosmological parameters are estimated by comparing the 
predicted signal m, to the data c/, as a function of scale Xj 
(which reduces to an angular scale a redhift z, or a mul- 
tipole {j respectively for lensing, supernovae and CMB data). 
We vary the parameters of the model, disposing them on a regu- 
larly spaced grid and evaluate, at each grid point, the likelihood 
function, 



£ = 



1 



(27!-)" | CI 1 / 2 



exp 



(29) 



Here C -1 is the data covariance matrix, including Poisson shot 
noise and cosmic variance. 

We focus on constraining the dark energy density and the 
parameter a of the quintessence potentials and restrict to a 
low-dimensional parameter space. Ideally, one would include 
at least nine cosmological parameters: the spatial curvature, the 
Hubble constant, the parameter a of the potential, the dark en- 
ergy and matter (both dark matter and baryonic) density param- 
eters, the reionization optical depth, the amplitude and spectral 
index of the initial power spectrum. In addition, one should 
include the three parameters accounting for the redshift distri- 
bution of sources, Eq. (I26> . 

Given the result of the analysis of CMB data 
JSpergel et al. 2003 1 we have assumed a spatially flat 
universe so that Q m o = 1 — Qqo- The amplitude of the initial 
power spectrum is fixed by the normalization on the CMB. 
We have also assumed that the reduced Hubble constant, the 
reionization optical depth, and the baryon energy density 
today are fixed to h = 0.72, T re i on = 0.17, and Qbo/z 2 = 0.024, 
respectively. 

We preliminary performed a Fisher matrix analysis on the 
parameters space (a, Qqo, n s ) in order to estimate approxi- 
mately the extent of the lcr region and decide the sampling 
steps of the grid for the computation of the likelihood. Figure[6] 
depicts the 68%, 95%, and 99% confidence levels (C.L.) af- 
ter for cosmic shear (top-hat variance) corresponding to a 
synthetic CFHTLS-wide like survey covering 170 deg 2 , with 
20 galaxies/arcmin 2 and an intrinsic ellipticity of 0.4. The top 
line refers to Ratra-Peebles models and the bottom line to 
SUGRA models. Interestingly, SUGRA models appear to be 
more constrained. Notice that the Fisher analysis holds only 
locally, around a fiducial model (marked by a cross) which is 
fixed at a = 6, Qqo = 0.73, and n s = 1.0. 

In the final likelihood analysis we allow to vary the cosmo- 
logical parameters in the following ranges 



Q Q o € [0.4,0.9], a e [0,25], n s e [0.9, 1.1]. 



(30) 



For the cosmic shear data, we also allow one of the source red- 
shift parameters to vary. We choose z s , always marginalized in 
the final analysis, over its 2cr interval; see Eq. J27> and OSt . 

Indeed, one can criticize these assumptions but these are 
sufficient i) to give us an idea of the parameter space avail- 
able for dark energy realized by quintessence and ii) to discuss 
to which extent weak lensing can improve the constraints on 
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Fig. 6. Fisher analysis of the top-hat shear variance of the cosmic shear on cosmological parameters (a, Oqo, n s ), for Ratra- 
Peebles (top line) and SUGRA (bottom line) models - contours at 68 and 95% C.L. We employ a synthetic realization of a 
CFHTLS-wide like survey, covering 170 deg 2 , with 20 gal/arcmin 2 and intrinsic ellipticity <r e = 0.4. The fiducial model is 
defined by (a, Qqo, n s ) - (6,0.73, 1) while (h, T re ion, Qboh 2 ) = (0.72,0.17,0.024) are kept fixed. The goal of this analysis is to 
help in choosing the range of the grid for the likelihood analysis. See § l3.4l for discussion. 



dark energy models. From a pragmatic point of view, we were 
limited by computational capacities and a more complete anal- 
ysis will follow. Let us stress that such an analysis requires an 
investigation of the potential degeneracies of the dark energy 
parameters with the standard cosmological parameters and in 
particular, one would need to quantify how allowing for dark 
energy changes the allowed range of variation of the other cos- 
mological parameters. This is left for further studies. 

3.5. Likelihood analysis: Joint cosmic shear data 

Using both real and synthetic cosmic shear data, we perform 
the likelihood analysis aiming to investigate to which extent the 
constraints on cosmological parameters, and in particular the 
quintessence ones, depend on the linear-to-non-linear map- 
ping and on the selection effects of the two-points statistics. 

In this section, we focus on the first issue by combining 
top-hat variance data of VIRMOS-Descart, CFHTLS-deep and 
CFHTLS-wide (22 deg 2 sub-sample) surveys; see Table^f™" 
numerical results of individual parameters constraints. Figure0 
depicts the results quoting the 68%, 95% and 99% confidence 
level contours. In particular, we compute the non-linear spectra 
by both the Peacock & Dodds (1996) and Smith et al. (2003) 
procedures, for Ratra-Peebles and SUGRA models as well, 
hence allowing for two kinds of comparisons. 

Firstly, we can compare Ratra-Peebles and SUGRA like- 
lihood contours, disregarding the non-linear mappings. There 
are two striking differences. One is the strong constraint on 
Qqo found in the SUGRA case. The reason for this is that, 
since the amplitude of the power spectrum is kept fixed, the 



(crg,fi m o) degeneracy implies a strong constraint on Q mu and 
consequently on Qqo through the flat universe prior. This ef- 
fect is stronger in SUGRA models than in Ratra-Peebles ones 
since the former approach a ACDM at low redshift. The other 
difference is the well defined degeneracy found in the (Qqo, a) 
plane for the Ratra-Peebles case. This feature will allow to put 
a stronger constraint on Ratra-Peebles' a than in SUGRA's a 
when combining with other data. 

Secondly and perhaps more interestingly, by comparing the 
non-linear mappings, for both Ratra-Peebles and SUGRA mod- 
els two classes of cosmological parameters come out: Those 
concerning the quintessence, a and Qqo, are found to be essen- 
tially independent of the non-linear mapping used. This sup- 
ports the claim of Simpson & Bridle (2005) that cosmic shear 
is sensitive to dark energy mostly through the background dy- 
namics. The second class involves the cosmological parameters 
accounting for the primordial universe, here n s only. In this case 
the corresponding likelihood contours strongly depend on the 
chosen mapping, regardless of the quintessence potential we 
used. Hence it seems not possible to constrain the primordial 
spectral index, at least jointly with quintessence parameters, 
until a stable formulation of the non-linear regime of structure 
formation will be available. 

3.6. Likelihood analysis: Synthetic data 

Using wide surveys, we can eventually investigate the effect 
of the non-linear regime of structures formation. For this pur- 
pose, we perform a likelihood analysis using a synthetic real- 
ization of the full CFHTLS-wide survey. This consists on syn- 



Carlo Schimd et al.: Tracking quintessence by cosmic shear 



11 




0.4 0.5 0.6 0.7 0.8 0.9 

n. 




1 


10 


1 


05 


c" 1 


00 





95 





90 


1 


10 


1 


05 


R* 1 


00 





95 





90 




0.4 0.5 0.6 0.7 0.8 0.9 



5 10 15 20 25 
a 



5 10 15 20 25 
a 




0.4 0.5 0.6 0.7 0.8 0.9 
fi 




1 


10 


1 


05 


d 1 


00 





95 





90 


1 


10 


1 


05 


R 1 


00 





95 





90 




5 10 15 20 25 
a 




0.4 0.5 0.6 0.7 0.8 0.9 

n. 




0.4 0.5 0.6 0.7 0.8 0.9 



1.10 
1.05 

R" 1.00 
0.95 



T 



0.90 I i 

0.4 0.5 0.6 0.7 0.8 0.9 
fi 





0.4 0.5 0.6 0.7 0.8 0. 
0. 



5 10 15 20 25 
a 



0.4 0.5 0.6 0.7 0.8 0. 



Fig. 7. Joint likelihood analysis of VIRMOS-Descart, CFHTLS-deep, and CFHTLS -wide/22 deg 2 top-hat variance data - con- 
tours at 68%, 95%, and 99% C.L. for the variables (Qoo,a), (a, n s ) and (£2qo, n s ). On the quintessence parameters sub-space, 
we have added the contour lines of (w p i vot , w a ) discussed in Fig. [3] assuming z p j vot = 0.5. The two upper lines are dedicated to 
Ratra-Peebles models [Eq. (|2ji] and the two lower to SUGRA models [Eq. @]. For each class of models, the non-linear spectrum 
has been computed using both the Peacock & Dodds (1996) procedure, first and third lines, and the halo model approach by 
Smith et al. (2003), second and fourth lines. See § !3.5l for discussion. 



thetic data vectors of top-hat and aperture mass variance and trix is computed using the analytical approximation derived in 

on a synthetic covariance matrix. The former are computed at |Schneider et al. (2002)| It depends on three main features of 

a fiducial model, which we take to be a ACDM with the cur- the survey; the effective area A, the effective galaxy number 

rent CFHTLS-wide redshift distribution. The covariance ma- density n ga i and the dispersion of the distribution of elliptici- 
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Fig. 8. CFHTLS-wide constraints on Ratra-Peebles models. From left to right, we present the likelihood analysis (68, 95, and 
99% C.L.) for the variables (Qqo, a), (Qqo, n s ) and (£2qo, a) and we have added the contour lines discussed in Fig. [5] assuming 
Zpivot = 0.5. The first line depicts the analysis of the actual dataset based on 22 deg 2 (W1+W3; see § I3.6> while the three other 
lines are based on the synthetic data for a field of 170 deg 2 . The second and third line show the top-hat shear variance and aperture 
mass variance, respectively. The fourth line describes the analysis of the same simulated data for the top-hat shear variance but 
using only angular scales larger than 20 arcmin to cut out the non-linear part of the matter power spectrum. In particular, we 
conclude from the left column that the parameters describing the quintessence sector (a and Oqo) are not affected by the choice 
of the statistics and are well estimated by the linear part of the power spectrum. Here we use the Peacock & Dodds (1996) 
mapping for the non-linear power spectrum. See § l3.6l for discussion. 



ties cr e . For these parameters we used the values that are ex- 
pected at the end of the CFHTLS-wide campaign, respectively: 
A = 170 deg 2 , n ga i = 20 gal/arcmin 2 , and cr e = 0.4. Notice 
we assume a larger ellipticity dispersion and a higher density 
of galaxies than those obtained in Hoekstra et al. (2005). They 



correspond to a different galaxy weighting scheme than the one 
used with current data. 

Figures [8] and [9] outline the likelihood analysis of Ratra- 
Peebles and SUGRA models, respectively - contours at 68%, 
95% and 99% confidence level. We show only the results 
achieved using the mapping by Peacock & Dodds (1996), those 
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Fig. 9. CFHTLS-wide constraints on SUGRA models. This figure is analogous to Fig. [8] We reach the same conclusion as for 
Fig-El We recover that the contours are almost independent on the value of a. See § B.6l for discussion. 



achieved using the mapping by Smith et al. (2003) being in 
agreement at a 10% level. 

Firstly, by comparing the results achieved using the top-hat 
variance data of the 22 deg 2 sub-sample (first line) with those of 
the synthetic 170 deg 2 field (second line), it is evident the gain 
achievable by the full survey. In particular, it is worth noticing 
that the pure quintessence parameters sub-space is less depen- 
dent on the survey area, while the constraints on the primordial 
spectral index strongly depend on it. Hence, the distinction of 



cosmological parameters in two classes introduced in the pre- 
vious section seems confirmed. 

The second and third lines show the likelihood contours 
for the top-hat and aperture mass variances, respectively. They 
are consistent, depicting compatible confidence level regions, 
for both Ratra-Peebles and SUGRA models. There are two 
main differences in the properties of aperture mass and top-hat 
variances. Measurements of aperture mass variance at differ- 
ent scales are less correlated than top-hat variance ones, since 
the former is a narrow filtered version of the shear power spec- 
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Fig. 10. Joint analysis of the quintessence parameters (£2qo, a, ) from CMB, Sn la and cosmic shear (top-hat variance), for Ratra- 
Peebles models (left panel) and SUGRA models (right panel). The likelihood analysis uses the "gold" set for Sn la (green 
contours), VIRMOS-Descart + CFTHLS-deep + CFTHLS-wide(22 deg 2 ) top-hat variance data for the cosmic shear (blue con- 
tours), and both combined (red contours). Contours correspond to 68, 95, and 99% C.L. According to WMAP-lyr measurements, 
the location of the first acoustic peak of the TT power spectrum of CMB, allowing for binning between neighboring multipoles, 
excludes regions of the parameter space (shadowed) nearly degenerate with Sn la constraints; see § !4.2l for details. On the right 
axes, we quote an indicative mass scale of the quintessence potential, weakly dependent on £2qo but Qqo is small; see §|2]for 
details. Finally, contours of w p i vot (solid) and w a (dotted) are superposed, setting z p i vot = 0. As for Ratra-Peebles models, super- 
novae data put strong constraints on a which are tighten by cosmic shear, while for SUGRA models all observables are fairly 
insensitive to the value of this parameter. The corresponding constraints are given on Tabled 



trum; and their effective range extend to only about 1/5 of the 
top-hat variance one. Hence they follow more accurately the 
shape of the power spectrum at their measured range. For a 
power spectrum featureless at most scales but these, the inde- 
pendent measurements render this statistic the most convenient 
to use. In this case, contours from aperture mass variance are 
expected to be smaller than top-hat variance ones. But in gen- 
eral, the loss of information regarding the behaviour of models 
at other scales diminishes the capability of distinguishing be- 
tween models. Consequently, the aperture mass variance con- 
tours are, in general, expected to be larger than top-hat variance 
ones. The difference is larger when analysing cosmic shear 
most important parameters. Since we are studying parameters 
to which cosmic shear is moderately sensitive to, the difference 
is not noticeable in our results of Fig.[S]and|5] but it is striking 
when using VIRMOS-Descart data to constrain the (<x 8 ,Q m ) 
plane (Ivan Waerbeke et al. 20 01 ). Furthermore, for our data 
with measurements up to less than 1 deg, the largest power 
spectrum scale probed by the aperture mass variance is the one 
probed by the top-hat variance at around 10 arcmin and its full 
range effectively lies on non-linear scales, rendering its infered 
parameters' constraints less reliable, due to non-linear model- 
ing uncertainties. For these two reasons, we will choose to use 
top-hat variance only, when producing the final results. 

By using the full wide survey, we can try to disentangle the 
effects of the non-linear regime of structures formation by cut- 
ting off the small angular scales from the final analysis. In such 
a way, we can better investigate the distinction of cosmologi- 
cal parameters in the two classes discussed above, probing if 
dark energy primarly hangs on the background dynamics. The 
plots on the bottom line of Fig. and [5] depict the analysis 
of the top-hat variance when taking into account only angu- 



lar scales larger than 20 arcmin, corresponding to wavevectors 
k ^ Ih Mpc~' atz S (z}/2 ~ 0.5, where therefore the ef- 
fects of the non-linear regime are sub-dominant. The spread 
in the likelihood contours is more relevant for plots involv- 
ing the primordial spectral index n s , while the quintessence 
parameters (a, Qqo) are not so much affected. Hence, one can 
study dark energy by cosmic shear using this technique even if 
not properly knowing how to deal with the non-linear regime. 
Obviously this conclusion has to be confirmed by a more com- 
plete study involving a larger parameter space, to account for 
other degeneracies. 

The first line of Fig. [8] and [9] show the current results of 
CFHTLS-wide, using Wl and W3 data from an effective sky 
coverage of 22 deg 2 . The corresponding marginalized con- 
straints on a and Oqo are a factor of 2 larger than the ones 
found with the synthetic data, depicted on the second line of 
the same figures. This is consistent with the reduction of the 
data error bars, which are proportional to o^n~^ A~ 1 * 2 , if cos- 
mic variance is not taken into account. We must caution that the 
gain in the parameters space cannot be estimated with precision 
by this simple argument, namely by assuming 



gain 



(31) 



where the several factors r are the ratios between the features 
of the two surveys. In fact, the gain in the data error bars does 
not translate linearly into a gain in the parameters space con- 
fidence levels. That happens only in the Fisher matrix approx- 
imation, and even there, only in the case of statistical uncor- 
rected parameters. Furthermore, this reasoning does not take 
into account the extra constraining power coming from mea- 
surements at larger scales, as discussed earlier on, or simply 
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Table 1. Results of separate and joint analysis of cosmic-shear and Sn la data on quintessence parameters (a, £2qo) at 68% (95%) 
confidence level. For Ratra-Peebles models, we quote only upper limits for a, the best-fit being always at a = 0. Only for 
indicative purposes, we quote 68% confidence level limits on (w p j vot , w a ) parameters computed at z p i vot = 0, setting n.c. when not 
constrained. Remind that in this class of models wq > -1. See § l4.1l for discussion. 





Ratra-Peebles 






SUGRA 




a 
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Q.QQ Wq 
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u '° -0.17 (-0.20) 
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n.c. 


n +12(+12) 
1J -5(-13) 
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n 7 .+0.09(+0.16) 
-O.lO(-O.ll) 
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U - -0.04(-0.04) 
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-0.13 


7 +7(+18) 

-K-2) 


n n . +0.03 (+0.03) 
-0.04 (-0.05) 


<, 


-0.84 


> 


-0.38 



coming from the fact of disposing of more degrees of freedom 
for the likelihood calculations. 

It is also worthwhile to notice that, in all cases, (from Fig.0 
to Fig. [5}, a ACDM model with a Harrison-Zel'dovich spec- 
trum [(a, n s ) = (0, 1)] is compatible with the data for Qqo ~ 0.7 
at 99% confidence level. 

4. Combining with other observables 

4.1. Snla 

We combine the cosmic shear data by VIRMOS-Descart and 
CFHTLS-deep and -wide (22 deg 2 sub-sample) surveys with 
the type la supernovae "gold" set by Tonry et al. (2003). In 
particular, we evaluate confidence intervals considering cosmic 
shear and Sn la both separately and jointly. Indeed, since the 
distance modulus depends only on the background dynamics, 
we can restrict to the quintessence cosmological parameters a 
and Qqo. 

Figure[Kj]depicts both the independent and combined anal- 
ysis for Ratra-Peebles and SUGRA models, using the top-hat 
variance data for the cosmic shear and computing the non- 
linear spectrum by the Peacock & Dodds (1996) mapping. 
The corresponding results are summarized in Table ^ these 
marginalized results for each parameter do not assume any 
prior knowledge of the other, apart from the flat priors implied 
by the range of the grid. Let us emphasize that the constraints 
obtained on (wq, w„) inferred from those on the parameters of 
the potential differ when we change the form of the potential. 
This confirms that constraints on the equation of state derived 
from a general parameterization have to be interpreted with 
care. 

Concerning Ratra-Peebles models, the weak lensing and 
Sn la contours are closer to mutual orthogonality, so nar- 
row joint constraints are expected. Data strongly favor a 
quintessence component close to a cosmological constant, the 
best-fit lying always at a = 0, hence in Table we prefer to 
present the constraints as an upper limit. As it is clear from 
Fig. [K)] Sn la are much more constraining for this parameter 
than weak lensing. However, even though Sn la alone reject a > 
1 at the 68% level, as was well known (Podariu & Ratra 2000 1, 
the information of weak lensing further narrows the interval. 

As far as SUGRA models are concerned, for a wide in- 
terval of a (approximately a > 5) both the luminosity dis- 
tance and the shear two-point correlations are (almost) inde- 
pendent of this parameter, as it is clear from Fig. [Hi] In fact, 



this conclusion was already reached studying the CMB temper- 
ature anisotropies (Brax et al. 2000> . arguing that the equation 
of state does not strongly depends on the slope of the poten- 
tial leaving both distances and linear growth factor almost un- 
changed. In terms of statistical significance, this means that the 
likelihood with respect to Sn la or weak lensing data alone is 
almost flat and so the best fit value has little meaning (see also 
ICaresia, Matarrese & M oscardini 2004 1. In Table [Owe list the 
results in the form of a parameter value, which we take to be the 
likelihood weighted average of a, plus or minus the necessary 
deltas to form the confidence intervals obtained. However, for 
a substantial interval of a, weak lensing and Sn la lead to dif- 
ferent, almost non-intersecting, ranges of Oqo: Sn la favoring a 
higher value of £2qo while weak lensing a lower one. Thus, even 
though both observables have limited sensitivity to constrain 
the parameter a, this fact allows to obtain a reasonable con- 
straint from the joint likelihood, as that interval of a is rejected 
by the data. The joint confidence interval are substantially re- 
duced, with a distinctive maximum of likelihood at a = 2. 
Notice, however, that this is an unstable situation. In fact, if 
cosmic shear or supernovae contours slightly change their ori- 
entation and size, due for example to a larger uncertainty on the 
redshift of sources, the joint contour will easily degradate. On 
the contrary, the extreme situation, with both contours vertical 
and parallel for all a, would imply the abandon of the subjacent 
quintessence model. It has to be stressed that every systematic 
effect on weak-lensing or supernovae data as well, relying for 
instance on the data analysis procedure, strongly affects the fi- 
nal result; therefore special care is necessary when combining 
several datasets. 



Provided that it is not evident which pivot redshift should 
be used when parameterizing these quintessence models by an 
equation of state of the form of Eq. (I12> . we superpose to the 
(Qqo, a) plane, in Fig. (1101 . the contours for w p j vo t and w a cor- 
responding to Zpi vot = 0. Let us stress that the final results in 
terms of the equation of state parameters are only indicative. 
Interestingly, notice from Fig. [3] that for Ratra-Peebles mod- 
els the estimation of w p i vot and w a should not change if using 



<.pivot 



0.5 instead of z 



pivot 



0. In fact, for this class of mod- 



els, noticeably different values on these parameters only appear 
when considering a high pivot redshift. On the contrary, for 
SUGRA models the choice of the pivot redshift would be rele- 
vant already at low redshift. See Table[2for specific constraints 
on wo and w a . 
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Fig. 11. Top line: Location of the first acoustic peak of the temperature CMB power spectrum, for Ratra-Peebles (left panel) 
and SUGRA (right panel) models. It is estimated using an analytical approximation lying on the solution of the background 
equations only (solid line), or by the full computation of Q s (dashed line). The central set of lines corresponds to the location 
of the first peak according to the best-fit of WMAP-lyr data, I - 220.1 ± 0.8 (Page et al. 2003 1, while the left and right ones 
correspond to the multipoles t — 201 and t = 240, defining the smallest and the largest multipoles contributing to the three points 
of binned data defining the first peak ( Hinsh aw et al. 2 003 1. Bottom line: Contour levels of cr% for Ratra-Peebles (left panel) and 
SUGRA (right panel) models with n s = 1, superposed to the cosmic shear data contours (VIRMOS-Descart + CFHTLS-deep + 
CFHTLS-wide/22deg 2 ). For indicative purposes, the star point marks the best-fit from the joint cosmic shear and SN la analysis, 
roughly sitting at cr g = 0.8. Remind that the spectrum was normalized on the CMB so this plots show the consistency with a 
normalization at z = 0. See § l4.2l for details. 



Recently, weak lensing data ( Semb oloni et al. 20051 
IHoekstra et a~ 2005 1 were used to put constraints on w under 
the assumption it is a constant parameter. In that particular 
case, one has two characteristic redshifts, z a and Zd e , defining 
the beginning of the acceleration phase [a(z a ) - 0], and of the 
domination of the dark energy [Q m (zde) = ^defoie)]- They are 
given by 

M_i_,-\3w_ 1 QdeO >3 W _ QdeO 

(1+Za) r— — — , (1+Zde) =7: — (32) 

1 + 3w i2 m0 I2 m0 

so that z a > Zde- When w becomes negative and large in ab- 
solute value, z a and za e tend to zero so that dark energy just 
starts almost today to dominate the universe and is redshifted 
in a way that it does not affect even low redshift observables. 
It follows that we expect the data to be insensitive to the value 
of w in that regime so that one can get only an upper limit on 
its value, as found by Semboloni et al. (2005) and Hoekstra et 
al. (2005). Indeed, such a situation cannot be achieved with a 
physical model as considered here because by construction it 



imposes that w > - 1 . It follows that our approach is a physi- 
cally motivated way of imposing a prior on w. Note also that 
Fig-fflshows that w = constant is not a good approximation of 
these two classes of models. At best the constraint on a con- 
stant w can be related to some redshift average of the equation 
of state. For these reasons, it is difficult to deduce a constraint 
on the physical models from a constraint on a constant w, even 
though the results of Tableware compatible with them. 

4.2. CMB 

CMB temperature anisotropics have been extensively 
used (see e.g. |Baccigalupi et al. 2002| IJassal et al. 20041 
Corasaniti et al. 20041 but several degeneracies amongst the 
cosmological parameters prevent to accurately constrain 
the cosmological parameters using CMB data only. Hence, 
concerning weak lensing, several studies already attempt to 
combine CMB and cosmic shear in order to constrain the 
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Fig. 12. Fisher analysis for Ratra-Peebles (left) and SUGRA (right) models, for cosmic shear top-hat variance (blue ellipses), 
Sn la "gold" set (green) and temperature CMB angular power spectrum by WMAP-lyr data (red) - contours at 68%, 95%, 
and 99% C.L. For the CMB, the noise matrix was computed using the public code by Verde et al. (2003) For cosmic shear, we 
consider both CFHTL-wide survey (A = 170 deg 2 , n ga i = 20 gal/arcmin 2 , cr e - 0.4) and a space based survey (A = 20, 000 deg 2 , 
«gai = 35 gal/arcmin 2 , <j e = 0.3), quoted as type II in the text. We consider only the quintessence parameters while the other 
parameters are kept fixed. See §|5]for discussion. 



cosmological parameters (see e.g. Contaldi, Hoekstra & Lewis 
2003; Ishak et al. 2004; Tereno et al. 2004). 

As for dark energy, it affects the CMB anisotropics angular 
power spectrum at least in two ways JBrax et al. 2000> . Firstly, 
the angular diameter distance is modified so that the peak struc- 
ture is shifted. In particular the location of the first acoustic 
peak, depending on the geometry of the universe, provides an 
estimate of the angular diameter distance to the last scattering 
surface. However, also pre-recombination effects can shift the 
peaks from their true geometrical locations dDoran et al. 2 000 
Kamionkowski & Buchalter 2000 1. Secondly, the time evolu- 
tion of the dark energy strongly affects the integrated Sachs- 
Wolfe effect. This effect is more relevant at low multipoles, 
modifying the amplitude of the spectrum, but it also leads to 
an additional shift of the Doppler peaks. 

Beside using CMB data to normalize the power spectra (see 
§ 13. 4> . we use them in still another way, by noticing that even 
without a statistical analysis of CMB data, we can strongly 
constrain cosmological parameters, and notably quintessence 
ones, simply by using the location of the first Doppler peak. 
Moreover, even without solving the perturbations equations 
and computing the TT spectrum, one can compute the acoustic 
scale just by solving the equation for the background evolu- 
tion. The location of acoustic peaks is then estimated allowing 
for the shifts induced by the dark energy by means of a fitting 
formulae ( |Doran & Lilley 2002) . In conclusion, the location of 
the first acoustic peak will be a function of the quintessence 
parameters (Qqo, a), with a negligible dependence on the pri- 
mordial spectral index n s due to the shift correction. Using 
this analytic approximation, we individuate, in Fig. [TO] the re- 
gion of the (Qqo, a) plane compatible with the location of the 
first acoustic peak of WMAP-lyr data (Hins haw et al. 2 003 
Page et al. 2003) including the bins' contributions to the three 
points defining the peak, namely 201 < £ < 240. Consistently 
with the normalization procedure we used, this result is in ex- 
cellent agreement with that deduced by the complete computa- 



tion of the TT spectrum, as shown in Fig.[^(top panels). Let 
us notice that the region of the parameter space (Qqo, a) com- 
patible with the position of the first acoustic peak is degenerate 
with the Sn la constraints; the reason being that, since the pre- 
recombination effects of quintessence on the definition of the 
acoustic horizon at last scattering are negligible, the acoustic 
scale eventually depends only on on the Hubble parameter like 
the luminosity distance, hence both are ultimately affected by 
quintessence approximatively in the same way. 

As mentioned before, we may compute, a posteriori, the 
cr 8 value for each CMB normalized model. The obtained cr 8 
iso-contours in the (Qqo, a) plane are shown in Fig. II lK bot- 
tom panels), setting n s = 1 . The <x 8 range that corresponds to 
the confidence levels found in this work, is in agreement with 
current cosmic shear <x 8 constraints; see Hoekstra et al. (2005) 
for a recent result and van Waerbeke & Mellier (2003) for a 
compilation of results. This shows the normalizations on the 
CMB and at z — are compatible. It is interesting to note that 
the directions of the Qqo - a degeneracy for lensing data are 
lines of constant <x 8 . Being the curvature of the universe kept 
fixed, a strong constraint in cr% implies a strong constraint in 
Qqo, through the well-known <x 8 - Q m o degeneracy, but only if 
Qqo and a are not much correlated. This is what we observe in 
the SUGRA case. In general, a contour in the (<x 8 , Q m ) plane 
will move as a function of a and a strong constraint in <x 8 does 
not necessarily imply a strong constraint in Qqo. This is what 
is obtained for the Ratra-Peebles case. 

The computational tools we have developed 
(see IRiazuelo & U zan 20001 ISchimd et al. 2005> allow us 
to compute, in the same framework, distance modulus, CMB 
anisotropies and weak lensing - cosmic shear effects. Hence, 
in principle there is no problem to combine Sn la and lensing 
data with CMB data. As already noticed, we have restricted 
here to a small parameter space. Indeed, a joint analysis with 
CMB data require a wider parameter space, possibly using the 
so-called "normal parameters" (see e.g. ISandvik et al. 20 04 1, 
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Fig. 13. Fisher analysis of weak lensing (top-hat variance) for Ratra-Peebles models, considering a CFHTLS-wide like survey 
and two space-based surveys layout (68, 95%, and 99% C.L.). Upper panels, analogous to the the likelihood analysis of Sec. [3] 
which also assumes all the other parameters fixed, compare a CFHTLS-wide like survey (wider ellipses) with a deep space-based 
survey with A = 1, 000 deg 2 , n ga \ = 50 gal/arcmin 2 , <r e - 0.3. Middle panels compare the deep space-based survey with a wider 
but shallow space-based survey with A = 20, 000 deg 2 , n ga i = 35 gal/arcmin 2 , cr e - 0.3. Bottom panels compare the CFHTLS- 
wide like survey with the type I space-based survey, but marginalizing over T re i on , left to vary. Compared with the first line, the 
last one shows the effect of a wider, more realistic parameters space. See §|5]for discussion. 



an option left for a future study. In such a case, we would 
integrate our pipeline with a Markov chain Monte Carlo code 
JTereno et al. 2 005 I developed for the likelihood analysis of 
cosmic shear. Furthermore, we will be able to include also the 
analysis of nucleosynthesis constraints, by means of a suitable 
code for quintessence models ( C oc et al. 2 006 1 which is able 
to deal with ordinary and extended quintessence models like 
the Boltzmann code and the lensing code used for this study. 

5. Weak lensing: prospects on future data 

Beyond CFHTLS, several next generation cosmic shear sur- 
veys are proposed in order to pin down w and the dark energy 
properties with exquisite details. It is interesting to explore the 
capabilities of such surveys to constrain the models discussed 
in this paper. 

A comprehensive study of all projects and observing strate- 
gies is however beyond the scope of this work, so we deliber- 
ately focus on rather simple concepts where both sky coverage 
and depth are increased, assuming systematics related to shape 



measurement can be discarded. In particular, we consider two 
possible layouts achievable by space-based missions. We will 
indicate by 

- type I, a deep survey that would cover about 1,000 deg 2 , 
one magnitude deeper than the CFHTLS-wide providing 
an effective galaxy number density around 50 gal/arcmin 2 ; 

- type II, a wider but shallow survey, covering 20,000 deg 2 
and yielding 35 gal/arcmin 2 . 

Furthermore, we assume an intrinsic ellipticity distribution of 
galaxies similar to the one observed with current cosmic shear 
surveys with HST. 

A Fisher analysis of Ratra-Peebles and SUGRA mod- 
els, restricted to the parameters (£2qo, a) keeping fixed the 
others, was performed around a fiducial model defined by 
(a, Q Q0 , n„ h, r reion , Q b a iy on()/i 2 ) = (6, 0.73, 1, 0.72, 0.17, 0.024). 
The results are presented in Fig. El Weak lensing surveys 
(blue ellipses) together with Sn la "gold" set (green ellipses) 
and CMB WMAP-lyr data (red ellipses) show that both 
quintessence parameters can in principle be determined with a 
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10% accuracy. The degeneracy with respect to Sn la and CMB 
would be almost totally broken. Is is worth noticing that, since 
all cosmological parameters but the quintessence ones are kept 
fixed, we have to take with care the astonishing gain with re- 
spect to a CFHTLS-wide like survey (larger ellipses) achiev- 
able by a space mission of type II (smaller ellipses). 

On scales larger than 10 degrees, the flat sky approxima- 
tion used in § B.ll does not hold anymore. Indeed, an angular 
distance of 15 degrees on a sphere has a 1% deviation from 
the same distance on a plane, and for larger scales spherical 
harmonics must be considered (Stebbins 1996 1. These effects 
are not taken into account in the results of Fig.^J where the 
gain observed between the two weak lensing ellipses essen- 
tially corresponds to Eq. i29\ . together with contributions from 
the fact that measurements from two surveys observing at dif- 
ferent scale ranges have different cosmic variances and dif- 
ferent degrees of freedom when fitting models to the data. In 
particular, a strategy of dividing the covered surveys areas in 
patches of 100 deg 2 is assumed. Hence, given the large ratio 
between both sky coverages, this is the dominant factor in the 
gain. 

For indicative purposes, and allowing for the effect of sys- 
tematic discussed in § 14.11 we restrict the Fisher analysis of 
weak lensing (top-hat variance) to Ratra-Peebles models, con- 
sidering a CFHTLS-wide like survey and both the space-based 
surveys of type I and type II; see Fig. ^] As in § 14.11 we 
evaluate the parameters space (or, £2qo, n s ). In the upper line, 
analogous to the likelihood analysis of Sec. [3] which also as- 
sumes all the other cosmological parameters fixed, we compare 
the CFHTLS-wide like survey (wider ellipses) with the space- 
based surveys of type I (smaller ellipses). The middle line is 
analogue, but comparing space-based surveys of type I (wider 
ellipses) with the space-based surveys of type II (smaller el- 
lipses); notice that the apparent rotation of the ellipses is simply 
due to a rescaling of the axes. As for quintessence parameters, 
a space-based survey of type II gives 99% C.L. contours ap- 
proximatively 3 times smaller than those achievable by a type I 
survey. Indeed more cosmological parameters have to be taken 
into account. As an example, allowing the reionization opti- 
cal depth T re ion to vary, the likelihood contours get strongly 
modified; see the third line of Fig. [D] where we compare a 
CFHTLS-wide like and type II space-based surveys, marginal- 
izing over Treion- They get larger and the degeneracies direc- 
tions are changed. The strong impact of the reionization optical 
depth in cosmic shear results comes from the degeneracy with 
between T re i on and the normalization of the spectrum. 

The great predictive power of space-based surveys will al- 
low to simultaneously constrain a large number of parameters 
with a good precision, even if not so high as the one showed 
in the analysis of a small number of parameters. However, it is 
important to notice that we have not taken into account the pos- 
sible use of tomography based on a decomposition of lensing 
data into several lensed/source planes, nor any use of higher or- 
der statistics than the top-hat shear variance. So, even if several 
issues have been neglected in the present study, our conclusions 
are likely not over-optimistic. 



6. Conclusions 

In this article, we have investigated the constraints set by 
weak lensing, supernovae and CMB data on two families of 
quintessence models. 

In such a situation where a physical model is fully speci- 
fied, we can treat both the background and perturbation evo- 
lution without any ambiguity, in particular when dealing with 
high redshift data. Such an approach is thus complementary, as 
discussed in Section 2, to those based on a parameterization of 
the dark energy sector, in particular when trying to infer con- 
straints on a physical model from those on the parameters of the 
equation of state. This also enables to get rid of the pivot red- 
shift problem when combining different data sets. From a more 
theoretical point of view, such models cannot lead to an equa- 
tion of state w < — 1, contrary to an arbitrary parameterization. 
It is thus interesting to determine whether there is a tension in 
the data when such a physical constraint is imposed, as would 
be concluded from various studies indicating that w < -1 is 
favored. 

To achieve this task, we have used a set of numerical tools 
that allow to compute background, CMB and lensing signatures 
of a large class of cosmological models including quintessence 
and some extensions such as scalar-tensor theories. We have fo- 
cused our analysis on three cosmological parameters, the index 
of the primordial power spectrum and two parameters describ- 
ing the quintessence models, along with an extra-parameter for 
the sources distribution (see § 13.41 . Although one can criticize 
such a small parameter space, it is sufficient to give an idea of 
the parameter space available for quintessence models and to 
discuss how weak lensing data can improve the constraint on 
dark energy. This choice was also driven by numerical limita- 
tions but our analysis will be extended to a larger set of param- 
eters in a near future. 

We have normalized the initial power spectrum to the CMB 
so that erg is now a prediction of the models and is not used for 
normalization. Note also that we do not require to specify an 
analytical form for the transfer function. Weak lensing predic- 
tions are also sensitive to the linear to non-linear mapping and 
we have discussed the effect of such mapping on the constraints 
with care. In particular, we have shown that, while the parame- 
ters of the primordial spectrum are sensitive, those of the dark 
energy sector remain robust. We have also tested the possibility 
to cut the weak lensing data sets (such as the CFHTLS-wide) 
in order to reduce the influence of the non-linear regime. 

This analysis is the first one using CFHTLS data to study 
the dark energy and illustrates the complementarity of these ob- 
servations with other data sets. To finish, we have also forecast 
how space-based wide field imagers will improve our knowl- 
edge of dark energy. In particular, we have considered two pos- 
sible strategies, the first deeper and the second wider but shal- 
lower. The latter turns out to be more suited to track dark en- 
ergy as far as cosmic shear is concerned. The constraints on the 
two classes of quintessence models considered in this article 
are shown on Fig. [H)] and Tabled They can be summarized as 
follows: For a flat universe and a quintessence inverse power 
law potential with slope a, we get a < 1 and Qqo = 0.75^ ^ 
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at 95% confidence level, whereas a = 2±' 8 , Q Q0 = 0.74+°°3 
when including supergravity corrections. 

In the future, we plan to improve this analysis by first com- 
paring it to a similar analysis based on a parameterization of 
the equation of state, by enlarging the parameter space, by ad- 
dressing more carefully the problem of the redshift distribution 
of galaxies and by shifting from a grid method to an MCMC 
method. 
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